Approximately self-similar critical collapse in 2+1 dimensions 
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Critical collapse of a self-gravitating scalar field in a (2+l)-dimensional spacetime with negative 
cosmological constant seems to be dominated by a continuously self-similar solution of the field 
equations without cosmological constant. However, previous studies of linear perturbations in this 
background were inconclusive. We extend the continuously self-similar solutions to solutions of the 
field equations with negative cosmological constant, and analyse their linear perturbations. The 
extended solutions are characterized by a continuous parameter. A suitable choice of this parameter 
seems to improve the agreement with the numerical results. We also study the dynamics of the 
apparent horizon in the extended background. 
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Following the numerical work of Choptuik and Prefo- 
rms on the critical collapse of scalar matter field in 2+1 
dimensions (see also Ref. 0), there has been debate 
about understanding and analytically reproducing their 
results l&yLtiaiul- Garfinkle found a one-parameter (n) 
family of continuously self-similar (CSS) solutions, and 
proposed that one of these is the critical solution for 
scalar field collapse in 2+1 dimensional AdS spacetime 
Numerical comparisons suggest that the critical 
value is n — 4. Subsequently, Garfinkle and Gundlach 
P| performed the linear perturbation analysis in this 
background and found that the only solution exhibiting 
a single growing mode is the n = 2 solution. Because 
of this discrepancy, they characterized their work as 
'inconclusive'. A weak point of their approach is that it 
neglects the negative cosmological constant A, although 
the latter is essential for the existence of black hole 
solutions in three dimensions (the BTZ black hole 
0). This was motivated in Ref. |j| by the following 
arguments: i) self-similarity requires A = 0; ii) close to 
the singularity, the A contribution to the full solution is 
negligible. Although these arguments seem reasonable, 
we expect the cosmological constant to play a crucial role 
in black hole formation [l|. Therefore, the inclusion of A 
in the above analysis may solve the above contradiction 
on the critical value of n. 



The Garfinkle family of CSS solutions is 
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where A = 2 2( - 1 ~ q ^ q q 2 and c 2 = l-l/2q. These solutions 
satisfy the three-dimensional Einstein equations 



Gab - Ag afc = V a 0V t 



2<?a 6 (V(/)) 2 



(2) 



with A = 0. The source term in Eq. is the stress- 
energy tensor of the minimally-coupled massless scalar 
field <p. (Note that the scalar field in Eq. J5J and Ref. Q 
differ by a factor — 1/^/Itt.) The Garfinkle CSS solutions 
are singular at u = v = 0. If q is a positive integer n, 
the initial region u > 0, v > can be extended across 
the surface v = 0, which plays the role of an apparent 
horizon, (q = n will be assumed below). 

We first extend Eq. Q to solutions of the Einstein equa- 
tions with A < and then consider the perturbation 
analysis in this background. Since the cosmological con- 
stant breaks the self-similarity, the appropriate variables 
are a scaling variable, for instance u, and a similarity 
variable, which we choose as y = (v/u) n . The metric 
coefficients and the scalar field are expanded in terms of 
the dimensionless combination Au 4 " 
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r = ^f^i B = r + Au bn F(y) + . . . , 
o = \ ln(2. 9w ,) = a + ku in G{y) + 
<f> = O + Ku in H{y) + . . . , 
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where ro, <7o and 4>o are t ne background contributions in 
Eq. QJ. At each order in the expansion the functions F, 
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G, and H satisfy a system of second-order coupled ordi- 
nary differential equations. We only consider the trunca- 
tion of the expansion (JHJ) to the first-order. The relevant 
equations are 
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plus two first-order contraints that reduce the moduli 
space of the initial conditions. Below we briefly discuss 
the solutions of Eqs. Q)-©. More details will be given 
in Ref. y| ■ The solution of Eq. regular at the center 
y = 1 (u — v) is 



where 



F(y) 



A 



y 5 f(y)dy + a(l-y 6 ) 
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The constant a can be set to zero by the gauge transfor- 
mation 
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can be obtained from Eq. (J5J. The two independent 
solutions of the homogeneous equation are 



Hi = 3 + 2y 2 + 3z/ 4 , 



ff 2 



(3 + 2j/ 2 + 3y 4 )ln 



1 + y 
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The regular solution of the inhomogeneous equation is 
H = C x Hi + C 2 H 2 , (11) 

where 

Ci = j XH 2 +c0, C 2 = - j" XHi , (12) 
X = ^[(1-^ + 2(2 + 3^] .(13) 

Finally G(y) is obtained by integration of Eq. JBJ with 
the boundary condition G(l) = —F'(l) — 0. This condi- 
tion follows from one of the constraints and implies the 
absence of conical singularities: 



ff* 1 " V,"ls=o = ^ a r tU r^\ y=0 
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The functions F 1 , G and H are shown to be analytic in 



z = y 1 ! 1 ' 



This first-order extension of the Garfinkle solutions is not 
uniquely defined, as we have found a one-parameter (/3) 
family of regular, analytic in z, solutions: 

(F, G, H) = {F, G, H) + c/3(0, Gp, Hp) , (15) 

where G and H are the (3 = solutions, and 
G,g = -c(l - y) 2 (3 + 2y + 3y 2 ), ify = H x . It can 
be shown that a new integration constant will appear at 
each successive order in the A-expansion. Thus, there 
is a manifold of exact solutions of Eqs. asymp- 
totic to the Garfinkle solutions near the singularity u = 0. 

Let us determine the effect of the first-order A-corrections 
on the location of the apparent horizon. From the defi- 
nition of the apparent horizon (Vr) 2 = 0, we obtain, to 
first-order in A, 



y 



2—l/n 



(l-A^(y)) =0, 



where 



4>(y) = 2G~QF+ ( 1+y ^ F' . 



y 



(16) 



(17) 



For A = 0, the apparent horizon is the past light cone 
y = of the singularity u = v = 0. For A < 0, the 
behavior of the functions F and G near y = is 



F(y) ~ F(0) + 



.4 



4(6n- 1) 



l/r. 



G(y) ~ G(0) + ../y 

4(5n — l)(6n — 1) 



(18) 



where F(0) and G(0) are determined numerically. On 
the apparent horizon, Eqs. (|18|l imply 
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Therefore, the apparent horizon recedes into the region 



y 



l/n 



< and becomes spacelike. This feature is 



essentially due to the term F'/y in Eq. (|17fl and does 
not depend on (3. 

The linear perturbation analysis in this background can 
be performed by expanding r, a and <j> as 

r = r + Au 6n F(y) + eu 2n - 2kn [f (y) + Au in h{y)] , 
a = a + Au* n G(y) + eu- 2kn [g (y)+Au An gi(y)} , 
<j> = 0o + Au An H(y) + eu- 2kn [h Q (y) + A U 4 "/n(y)](20) 

where e is a small parameter that controls the strength of 
the perturbation, and we have truncated the expansion 
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to first-order in Au 4n . The growing modes are given by 
Re(k) > 0. The critical solutions have a single growing 
mode 0. 

The analysis of the zeroth-order perturbations fo, go and 
ho was carried out in Ref. £|. Here, we only recall the 
main points of this analysis. The regular solution of the 
differential equation for f (/o(l) = 0), 



is 



/o = Cl (l-2/ 2 - 2fc ). 



(21) 



(22) 



This solution is pure gauge, i.e. it can be generated from 
the unperturbed solution ro(u, v) — (u 2n — v 2n )/2 by the 
coordinate transformation 



V n 



-2kn 



(23) 



In the gauge fo = 0, the scalar field perturbation 
u~ n ho(y) solves the massless Klein-Gordon equation 
for the A = background spacetime. The solution, 
in terms of hypergeometric functions, depends on two 
integration constants. The first one is fixed by the 
regularity condition, i.e. the absence of logarithmic 
divergence for y = 1. The second integration constant 
is fixed by the condition of smoothness on the null line 
y = 0, i.e. the analyticity (in at least one gauge) of ho 
as a function of y 1 ^' 1 — v/u. A necessary condition is 
2kn = m, where m is a positive integer. For m < n this 
condition is also sufficient. For n < m < 2n one can find 
a gauge, i.e. a value of ci, such that ho is analytic. By 
contrast, for m = n and m > 2n there is no gauge in 
which ho is analytic. The second-order equation for go 
shows that go is generically divergent on the null line 
y = 0. However, for 1 < m < n there is a gauge in which 
go is analytic, and for the value m = 2n — 1 > n (n > 1) 
go and ho are analytic in the same gauge. Therefore, 
regularity at the origin and analyticity in y 1 /™ require 
k = m/2n and either i) 1 < m < n or ii) m = 2n — 1 
(n > 1). It can be easily seen that only the solution 
with n = 2 has a single unstable mode, namely m = 3 
(k = 3/4). 

The only debatable question in this analysis is whether 
the requirement that the non-scalar quantity go is ana- 
lytic at y = might not be too strong. In principle, it 
should be enough to demand that the perturbation of a 
scalar quantity, such as the Ricci scalar, is analytical at 
y = 0. At zeroth-order in A, the Ricci scalar is 



In the gauge c\ = 0, 50 an d h' diverge for y — > as 
Therefore, po diverges. However, the zeroth- 
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-m/n 



order Ricci scalar Rq goes to zero as y 1 l ' n and the per- 
turbation 



5R oc Ropo 
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-(m+1) / n 



(26) 



remains finite at y = for m < n, including m = 1. If 
the mode m = 1 (k = l/2n) were allowed, none of the 
Garfinklc solutions would be critical: for n = 1 there 
would be no growing mode, for n = 2 there would be 
two growing modes, m — 1 and m = 3, etc. However, 
as we now show, the extra modes with m = 1 do not 
survive the first-order extension in A. 

The first-order perturbation fx solves the inhomogeneous 
differential equation 

-yf»+(5-2k)fi = A^ ( l +y) ^^[ /o + ( l_ y 2 )5o ]. 

(27) 

In the gauge fo = 0, go ~ y m ^ n implies 

fi ~ y~^~ {m > 1) or hi V {m = 1) (28) 

for y — > 0. If m > 1, the divergence of fi can be 
gauged away by the zeroth-order gauge transformation 
(|23|l . The logarithmic divergence of the first-order 
contribution to the m = 1 mode cannot be gauged 
away; this mode is never analytic at y = 0. A detailed 
analysis of the first-order perturbations will be presented 
elsewhere |8j . Here, let us just note that the analytic and 
numerical integrations of the first-order perturbations 
indicate that all the modes found in the analysis of 
the zeroth-order perturbations satisfy the boundary 
conditions of regularity and analyticity at first-order for 
any value of j3. This shows that the analysis of Garfinklc 
and Gundlach is robust, i.e. it survives extension to 
first-order in A. 

Now let us discuss the effect of the extension on the be- 
havior of the apparent horizon for the perturbed critical 
solution (n = 2, k = 3/4). The apparent horizon satisfies 
the equation 



Y^yTj f 1 - Au V _ eu ~ 3 X - eAu5 vj = , (29) 



where ip is defined in Eq. (|17fl and 



fo- 



(30) 



R = Ro(u, y )(l-2eu- 2kn p Q (y)), 



where 



Po = go 



i + y 

4c 



[(1 - y)h' - 2kh ] 



(24) 



(25) 



(The exact form of 77 is inessential for the following dis- 
cussion and will be given in Ref. 0.) For y — > 0, x is 
dominated by the last term in Eq. (|30|) : 



.£1^-3/2 

2 y 



(31) 



FIG. 1: The scalar field as a function of R for the n = 2 CSS 
solution (solid), the n = 2 extended CSS solution (dashed), 
and the n — 4 CSS solution (dotted). The n = 2 extended 
CSS solution is computed for /3 = 1 and To = 1.9. The 
scalar field has been rescaled by a factor — l/y/4-ir to facilitate 
comparison with the results of Ref. 



The zeroth-order approximation of Eq. (|29|l with A = 
is 



eci 
' 2 



(32) 



The apparent horizon is null, and exists for both signs of 
e. (The singularity u = v = is hidden by the apparent 
horizon only for eci < 0.) The situation changes dra- 
matically when we take into account the first-order con- 
tributions. Neglecting the term rj, we see from Eq. I|19f) 
and Eq. I|31|l that near y = the shape of the apparent 
horizon is determined by a balance between the A and e 
contributions. The leading behavior is 



u 



V 3 



u 



m 1 " 



(33) 



The apparent horizon, which exists only for eci < 0, is 
spacelike for small positive y and becomes null (u = uq) 
for y — 0. The numerical solution of Eq. I|29l) shows that 
on the apparent horizon u is everywhere bounded by 
wo. This confirms a posteriori that the r\ contribution 
to Eq. I|29(l can be neglected for small e. The existence 
and the shape of the apparent horizon, which hides the 
singularity u — v = 0, do not depend on the parameter (3. 

Finally, we present some evidence that the O(A) correc- 
tions improve the agreement with the numerical simu- 
lations of near-critical collapse. Following Ref. we 
introduce the coordinates (T, R) 

1 - ii 2 

T=-2n\nu, R = U - 2n r = — + Ae~ 2T F(y) . 

(34) 

The expression of the extended Garfinkle scalar field at 



FIG. 2: The mass aspect M as a function of R for the 
n — 2 CSS solution (solid), the n = 2 extended CSS solu- 
tion (dashed), and the n — 4 CSS solution (dotted). 



some fixed T is 



<l>n(y,To) 



-2c In 



(i + y) 



Ae- 2T "(H(y)-H(l)), 

where has been shifted by a constant to make it vanish 
at y = 1. In Ref. Garfinkle shows that the nonex- 
tended solution with n = 4 agrees with the numerical 
critical solution of Ref. Q] at an intermediate time To ~ 9. 
For such a large To , the extended <j> n (|35|l reduces to that 
of the CSS solution. However, the calibration of the nu- 
merical To involves some ambiguity. In Ref. 0, T is 
defined by T = — ln£ c , where t c = (T = +oo) at the 
accumulation point (the singularity). Even a tiny er- 
ror in the determination of this zero from near-critical 
simulations will translate into a large error on the cor- 
responding value of To. Therefore, the latter has to be 
considered as an unknown parameter. A second unknown 
parameter in Eq. (|3~5|l is (3. T and (5 can be set by 
comparing 02(0, To) for the critical solution n = 2 with 
04(0, oo). From the numerical solution of Eqs. Q}-©, we 
find H(Q) - 0.016 < c(3H p (0) = 3c/3, provided that j3 
is not too small. So 02(0, T)) depends only on the prod- 
uct (3e~ 2T ° . By comparing the latter with the numerical 
value of 4 (O,oo), we obtain /3e" 2T ° - 0.022. For the 
(arbitrary) choice (3 = 1 this gives Tq ~ 1.9. In Fig. 1 
we plot in terms of R the CSS solution with n = 2, the 
extended solution for n = 2 and the CSS solution with 
n = 4. The O(A) corrections seem to improve the agree- 
ment of the n = 2 critical solution with the numerical 
results. This conclusion is strengthened by Fig. 2, where 
we plot the mass aspect in terms of R 



M n (y,T Q ) = -Ar 2 + Ae^r, u r, „ 
\l-y 2 ) 2 ( Ay 



4/y 



- Au 4 ™ 



(i + y) 2 

-l 

ip(y) 



(36) 
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for the same values of (3 and To- However, it is clear 
that the first-order extended n — 2 solution agrees with 
the numerics only over a small range of To, as opposed 
to the n = 4 CSS solutions, which agree over a large 
range of intermediate times [3|]. This suggests that the 
question of the agreement between the analytical and 
numerical critical solutions is still an open problem. 

To conclude, our analysis shows that in the near-critical 
regime the shape of the apparent horizon is determined 
by a balance betwen the A and e contributions. This is 
evidence that the cosmological constant plays a role in 
black hole formation. We have also shown that the ap- 



parent contradiction between the results of Ref. 1] and 
Ref. can partly be solved by including 0(A) terms. 
Another result of our analysis is that, at this order, there 
seems to be a one-parameter family of critical solutions, 
rather than a single critical solution. This parameter 
is not connected with gauge transformations (the gauge 
parameter is a, which has been set to zero). Rather, as 
will be discussed in more detail in Ref. , the first-order 
terms linear in {3 in Eq. I|15|) . which solve the homoge- 
neous equations can be reinterpreted as zeroth- 
order k = — 2 perturbations. It follows that the extended 
n = 2 critical solution is unique modulo the addition of 
a decaying perturbation. 
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